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\Q [ Abstract 

f— - n , Let P, Q C R 2 be two n-point multisets and > b be a set of A inequalities on x and y, where 

A e R Ax2 , r = [y], and 6 e R A . Define the constrained Minkowski sum {P®Q)Ar>b as the multiset 

Q ! {(p + <z)|p e P,q e Q, A(p + q) > b}. Given P, Q, > 6, an objective function / : R 2 -> R, and 

a positive integer fc, the Minkowski Sum Selection problem is to find the k th largest objective 
value among all objective values of points in (P®Q)Ar>b- Given P, Q, Ar > b, an objective function 
/ : R 2 — > R, and a real number S, the Minkowski Sum Finding problem is to find a point (x*,y*) 

,_l | in (P © Q)Ar>b such that \f(x*,y*) - S\ is minimized. For the Minkowski Sum Selection 

J> ■ problem with linear objective functions, we obtain the following results: (1) optimal 0(n log n) 

t-H I time algorithms for A = 1; (2) 0(nlog 2 n) time deterministic algorithms and expected Oinlogn) 

time randomized algorithms for any fixed A > 1. For the Minkowski Sum Finding problem with 
linear objective functions or objective functions of the form f(x,y) — jj^, we construct optimal 
• ■ O(nlogn) time algorithms for any fixed A > 1. As a byproduct, we obtain improved algorithms 

for the Length-Constrained Sum Selection problem and the Density Finding problem. 

■ Keywords. Bioinformatics, Sequence analysis, Minkowski sum. 

o : 

■ 1 Introduction 

Let P,QCM? be two n-point multisets and Ar > b be a set of A inequalities on x and y, where 
A G M Ax2 , r = [y], and b G R A . Define the constrained Minkowski sum (P © Q)Ar>b as the multiset 
{(p + q)\pe P,qeQ,A(p + q) > b}. 

In the Minkowski Sum Optimization problem, we are given P, Q, Ar > b, and an objective 
function / : R 2 — > R. The goal is to find the maximum objective value among all objective values 
of points in (P © Q)Ar>b- A function / : D C IR 2 — > R is said to be quasiconvex if and only if for 
all points t>i,f2 G -D and all 7 G [0, 1], one has 7(7 • V\ + (1 — 7) • ^2) < max(/(t>i), /(t^))- Bernholt 
ei aZ. [5j studied the Minkowski Sum Optimization problem for quasiconvex objective functions and 
showed that their results have applications to many optimization problems arising in computational 
biology PEl[3[Tni[I3[IS[T3[2T]. In this paper, two variations of the Minkowski Sum Optimization 
problem are studied: the Minkowski Sum Selection problem and the Minkowski Sum Finding 
problem. 

In the Minkowski Sum Selection problem, we are given P, Q, Ar > 6, an objective function / : 
R 2 — ► R, and a positive integer k. The goal is to find the k th largest objective value among all objective 
values of points in (P © Q)Ar>b- The Minkowski Sum Optimization problem is equivalent to the 
Minkowski Sum Selection problem with k = 1. A variety of selection problems, including the Sum 
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Selection problem [3J[T3], the Length-Constrained Sum Selection problem [J2], and the Slope 
Selection problem [18], are linear-time reducible to the Minkowski Sum Selection problem 
with a linear objective function or an objective function of the form f(x,y) = It is desirable that 
relevant selection problems from diverse fields are integrated into a single one, so we don't have to 
consider them separately. Next, let us look at the use of the Minkowski Sum Selection problem 
in practice. As mentioned above, the Minkowski Sum Optimization problem finds applications to 
many optimization problems arising in computational biology [D El [101 HI EE 13 [21] . In these 
optimization problems, the objective functions are chosen such that feasible solutions with higher 
objective values are "more likely" to be biologically meaningful. However, it is not guaranteed that 
the best feasible solution always satisfies the needs of biologists. If the best feasible solution does not 
interest biologists or does not provide enough information, we still have to find the second best feasible 
solution, the third best feasible solution and so on until a satisfying feasible solution is finally found. 
As a result, it is desirable to know how to dig out extra good feasible solutions in case that the best 
feasible solution is not sufficient. 

In the Minkowski Sum Finding problem, we are given P, Q, Ar > b, an objective function 
/ : M 2 —> R, and a real number 5. The goal is to find a point (x*,y*) in (P © Q)Ar>b such that 
\f(x*,y*) — 5\ is minimized. This problem originates from the study of the Density Finding problem 
proposed by Lee et al. [12] . The Density Finding problem can be regarded as a specialization of 
the Minkowski Sum Finding problem with objective function f(x,y) = ^ and find applications in 
recognizing promoters in DNA sequences |11[ 120] . In these applications, the goal is not to find the 
feasible solution with the highest objective value. Instead, feasible solutions with objective values 
close to some specific number, say 5, are thought to be more biologically meaningful and preferred. 

The main results obtained in this paper are as follows. 

- The Minkowski Sum Selection problem with one constraint and a linear objective function 
can be solved in optimal O(nlogn) time. 

- The Minkowski Sum Selection problem with two constraints and a linear objective function 
can be solved in 0(n log 2 n) time by a deterministic algorithm and expected 0(n log n) time by 
a randomized algorithm. 

- For any fixed A > 2, the Minkowski Sum Selection problem with A constraints and a linear 
objective function is shown to be asymptotically equivalent to the Minkowski Sum Selection 
problem with two constraints and a linear objective function. 

- The Minkowski Sum Finding problem with any fixed number of constraints can be solved in 
optimal O(ralogn) time if the objective function f(x,y) is linear or of the form — . 

As a byproduct, we obtain improved algorithms for the Length-Constrained Sum Selection 
problem [14] and the Density Finding problem [12] . Recently, Lin and Lee [J3] proposed an expected 
0(n log(u — 1 + l))-time randomized algorithm for the Length-Constrained Sum Selection prob- 
lem, where n is the size of the input instance and !,u6N are two given parameters with 1 < I < u < n. 
In this paper, we obtain a worst-case 0{n log(u — I + l))-time deterministic algorithm for the Length- 
Constrained Sum Selection problem (see Appendix A). Lee, Lin, and Lu [12] showed the Density 
Finding problem has a lower bound of Q(nlogn) and proposed an 0(nlog 2 m)-time algorithm for 
it, where n is the size of the input instance and m is a parameter whose value may be as large as n. 
In this paper, we give an optimal 0{n log n)-time algorithm for the Density Finding problem (see 
Appendix B). 
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2 Preliminaries 



In the following, we review some definitions and theorems. For more details, readers can refer to [5l[8]. 
A matrix X £ M nxm is said to be sorted if the values of each row and each column are in nondecreasing 
order. Frederickson and Johnson [8] gave some results about the selection problem and the ranking 
problem in a collection of sorted matrices. From the results of Frederickson and Johnson [8], we have 
the following theorems. 

Theorem 1: The selection problem in a collection of sorted matrices is given a rank k and a collection 
of sorted matrices {Xi, . . . , X^} in which Xj has dimensions rij x rrij, rij > rrij, to find the k th largest 
element among all elements of sorted matrices in {X\, . . . ,Xjy}. This problem is able to be solved in 
0(Y^f = i rrij log (2nj/nij)) time. 

Theorem 2 : The ranking problem in a collection of sorted matrices is given an element and a 
collection of sorted matrices {Xi, . . . , Xjy} in which Xj has dimensions rij x rrij, rij > rrij, to find the 
rank of the given element among all elements of sorted matrices in {X\, . . . ,Xn}. This problem can 
be solved in 0(J2jL\ rrij log(2n j /rrij)) time. 

By the recent works of Bernholt et al. [5], we have the following theorems. 

Theorem 3: Given a set of A linear inequalities Ar > b and two n-point multisets P, Q CM 2 , one 
can compute the vertices of the convex hull of (P © Q)Ar>b in 0(Alog A + A • nlogra) time. 

Theorem 4: The problem of maximizing a quasiconvex objective function / over the constrained 
Minkowski sum (P © Q)Ar>b requires fi(nlogn) time in the algebraic decision tree model even if / is 
a linear function and Ar > b consists of only one constraint. 

3 Minkowski Sum Selection with One Constraint 

In this section we study the Minkowski Sum Selection problem and give an optimal 0{n log n) 
time algorithm for the case where only one linear constraint is given and the objective function is also 
linear. 

3.1 Input Transformation 

Given P = {(2:1,1, yi,i), (x 1>n , yi, n )}, Q = {(^2,1,2/2,1), (x 2 ,n,y2,n)}, a positive integer k, one 
constraint L: ax + by > c, and a linear objective function f(x,y) = dx + ey, where a, b, c, d, and e 
are all constants, we perform the following transformation. 

1. Change the content of P and Q to {(0x1,1 + 62/1,1, <&Gi,i + ej/i,i), . . . , (axi, n + 6yi, n , da;i, n + eyi, n )}, 
and {(022,1 + 6j/ 2 ,l, dx 2 ,i + ey 2 ,i), ■■■ , (ax 2 , n + by 2 , n , dx 2)Tl + ey 2>n )}, respectively. 

2. Change the constraint from ax + by > c to x > c. 

3. Change the objective function from dx + ey to y. 

This transformation can be done in 0(n) time and the answer remains the same. Hence from now 
on, our goal becomes to find the k th largest y-coordinate on the constrained Minkowski sum of P and 
Q subject to the constraint L : x > c. 
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3.2 Algorithm 



For ease of exposition, we assume that no two points in P and Q have the same x-coordinate and n is a 
power of two. The algorithm proceeds as follows. First, we sort P and Q into P x and Q x [P y and Q y , 
respectively) in nondecreasing order of x-coordinates (y-coordinates, respectively) in 0(n log n) time. 
Next, we use a divide-and-conquer approach to store the y-coordinates of (P © Q) x > c as a collection 
of sorted matrices and then apply Theorem Q] to select the k th largest element from the elements of 
these sorted matrices. 

Now we explain how to store the y-coordinates of (P © Q) x >c as a collection of sorted matrices. 
Let P x = ((xi,yi),...,(x n ,y n )), Q x = ((§1, yi), . . . , (x n , y„)), P y = ((x' 1 ,y' 1 ), . . . , (x' n ,y' n )), and Q y = 
((x'l, y'i), • • • , (x' n ,y' n )). We then divide P x into two halves of equal size: A = ((xi, yi), ■ ■ ■ , (x n / 2 ,y n / 2 )) 
and B = ((x n / 2+1 , y n / 2 +i), ■■■ , (x n , Vn))- Find a point (x t , y t ) of Q x such that x n / 2 +x t < c and t is max- 
imized. Then divide Q x into two halves: C = ((xi,yi),... ,(x t ,y t )) and D = ((x t +l, Vt+i), ■ ■ ■ ,(x n ,y n )). 
The set (P © Q) x > c is the union of (A © C) x > c , (A © D) x > c , (B © C) x > c , and (B © D) x > c . Because x t 
is the largest x-coordinate among all x-coordinates of points in Q x such that x n / 2 + xt < c, we know 
that all points in A © C cannot satisfy the constraint x > c. Hence, we only need to consider points 
in A © D, B © C, and B ® D. Because P x and Q x are in nondecreasing order of x-coordinates, it is 
guaranteed that all points in B © D satisfy the constraint L, i.e., B © D = (B ® D) x > c . Construct 
in linear time row-vector = (ri,r 2 , ■ ■ ■ ,r n / 2 ) which is the y-coordinates in the subsequence of P y re- 
sulting from removing points with x-coordinates no greater than x n / 2 from P y . Construct in linear 
time column_vector = (c\,c 2 , . . . ,Cn-t) which is the y-coordinates in the subsequence of Q y resulting 
from removing points with x-coordinates no greater than xt from Q y . Note row-vector is the same 
as the result of sorting B into nondecreasing order of y-coordinates, and column-vector is the same 
as the result of sorting D into nondecreasing order of y-coordinates. Thus, we have {y : (x, y) G 
(B © D) x > c } = {y : (x, y) G B © D} = {ri + Cj : 1 < i < \rowjvector\, 1 < j < \column_vector\}. 
Therefore, we can store the y-coordinates of B © D = {B © D) x > c as a sorted matrix X of dimensions 
\row-vector\ x \columri-vector\ where the (i,j)-th element of X is rj + Cj. Note that it is not necessary 
to explicitly construct the sorted matrix X, which needs Q(n 2 ) time. Because the (i,j)-th element of 
X can be obtained by summing up and Cj, we only need to keep row -vector and column-vector . The 
rest is to construct the sorted matrices for the y-coordinates of points in (A © D) x > c and (B © C) x > c . 
It is accomplished by applying the above approach recursively. The pseudocode is shown in Figure CD 
and Figure [5J We now analyze the time complexity. 

Lemma 1: Given a matrix X G ]R 5r{xOT , we define the side length of X be 91 + SOT. Letting T(n', ml) 
be the running time of ConstructMatrices(P x , Q x , P y , Q y , L), where n' = \P X \ = \P y \ and m! = \Q X \ = 
\Qy\i we have T(n',m') = 0((n' + m')log(n' + 1)). Similarly, letting M(n',m') be the sum of the 
side lengths of all sorted matrices created by running ConstructMatrices(P r , Q x , P y , Q y , L), we have 
M(ri, to') = 0((ri + to') log(n' + 1)). 

Proof: It suffices to prove that T(n', to') = 0((n'+m') log(n'+l)). By Algorithm ConstructMatrices in 
Figured! we have T(n', to') < maxo<i< m '{c'(n'+TO')-|-T(n'/2, i)+T(n'/2, to'— i)} for some constant d . 
Then by induction on n', it is easy to prove that T(n', to') is 0((n' + to') log(n' + 1)). □ 

Theorem 5: Given two n-point multisets PCK 2 and Q C M 2 , a positive integer k, a linear constraint 
L, and a linear objective function / : M 2 — > R, Algorithm Selectioni finds the k th largest objective 
value among all objective values of points in (P © Q)i in 0(n log n) time. Hence, by Theorem [H 
Algorithm Selectioni is optimal. 
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Algorithm ConstructMatrices(Pr, Q x , P y , Q y , L: x > c) 

Input: P x and P y are the results of sorting the multiset P C R 2 in nondecreasing 
order of x-coordinates and y-coordinates, respectively. Q x and Q y are the results 
of sorting the multiset Q C M 2 in nondecreasing order of x-coordinates and 
^-coordinates, respectively. A linear constraint L: x > c. 

Output: The y-coordinates of points in (P © Q) x >c as a collection of sorted matrices. 

1 n' <- |Pa;|; m' <- \Q X \. 

2 Let P x = ((x 1 ,y 1 ),...,(x n >,y n >)), Q x = ((xj, yi), . . . , (x m >, y m ')), 
p y = {{x'i,y[),... ,(x' n ,,y' n ,)), and Q y = ((x' 1; y[), . . . , \x' m ,,y' m ,)). 



3 xo < oo. 

4 if n' < or m' < then 

5 return 

6 if n' = 1 or m! = 1 then 

7 Scan points in Pj, © to find all points satisfying L and construct the sorted 
matrix for y-coordinates of these points. 

8 return the above sorted matrix. 

9 for t <— m! down to do 

10 if x n i j2 + xt < c then 

11 row-vector <— subsequence of P y being removed points with 
x-coordinates < x n / 

12 column_vector <— subsequence of Q y being removed points with 
x-coordinates < Xf 

13 A x «- P x [l,n'/2]; B x «- P x [n'/2 + 1, n']; C s «- Q«[l,*]; D a «- + l,m']. 

14 A y <— subsequence of P y being removed points with x-coordinates > x n // 2 - 

15 B y <— subsequence of P y being removed points with x-coordinates < x n // 2 - 

16 C y <— subsequence of Q y being removed points with x-coordinates > Xf 

17 D y <— subsequence of Q y being removed points with x-coordinates < x%. 

18 ConstructMatrices(A x , D x , A y , D y , L: x > c). 

19 ConstructMatrices(Pa;, C x , B y , C y , L: x > c). 

20 return row-vector and column_vector. 



Figure 1: The subroutine for the Minkowski Sum Selection problem with one linear constraint 
and a linear objective function. 

Proof: Let S = {X\, . . . , Xjy} be the sorted matrices produced at Step 4 in Algorithm Selectioni. Let 
Xj, 1 < j < N, be of dimensions rij x rrij where rij > rrij. By Lemma[TJ we have OQ^- ,(mj + = 
O(nlogn). By Theorem [IJ the time required to find the k th largest element from the elements of 
matrices in S is 0(Xa=i m « log(2r^/mj)). Since 

jv ^271^^ 

^mjlog(2ni/mj) < y^mj — 1 = ^2n; < ^2(m; + rij) = 0(n log n), 

i=l i=l 1 i=l i=l 

the time for selecting the k th largest element from elements of matrices in S is 0(n log n). Combining 
this with the time for the input transformation, sorting, and executing ConstructMatrices(P x , Q x , P y , Q y , L), 
we conclude that the total running time is 0(n log n). □ 
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Algorithm Selection! (P, Q, L, /, k) 

Input: Two multisets PCM 2 and Q C R 2 ; a linear constraint L; a linear objective 

function / : R 2 — > R; a positive integer fe. 
Output: The largest value among all objective values of points in (P © Q)l- 

1 Perform the input transformation described in Section 3.1. 

2 Sort P and Q into P x and Q x , respectively, in nondecreasing order of x-coordinates. 

3 Sort P and Q into P y and Qy, respectively, in nondecreasing order of y-coordinates. 

4 S <— ConstructMatrices(P x , Qa;, Py, Qy, L). 

5 return the fc*' 1 largest element among the elements of sorted matrices in S. 



Figure 2: The main procedure for the Minkowski Sum Selection problem with one linear constraint 
and a linear objective function. 

Using similar techniques, the following problem can also be solved in O(nlogn) time. Given two 
n-point multisets P C R 2 and Q C R 2 , a linear constraint L, a linear objective function / : R 2 — > R, 
and a real number t, the problem is to find the rank of t among all objective values of points in 
(P © Q)l, where the rank of t is equal to the number of elements in {y\(x,y) 6 (P © Q)l,U > t} 
plus one. The pseudocode is given in Figure El Note that in Algorithm Selectioni and Algorithm 
Rankingi, we assume the input constraint is of the form ax + by > c. After slight modifications, we 
can also cope with constraints of the form ax + by > c. To avoid redundancy, we omit the details 
here. For ease of exposition, we assume that Algorithm Selectioni and Algorithm Rankingi are also 
capable of coping with constraints of the form ax + by > cm. the following sections. 



Algorithm Rankingi (P, Q, L, /, t) 

Input: Two multisets P C R 2 and Q C R 2 ; a linear constraint L; a linear objective 

function / : R 2 — > R; a real number t. 
Output: The rank of t among the objective values of points in (P © Q)l- 

1 Perform the input transformation in Section 3.1. 

2 Sort P and Q into P x and Q x , respectively, in nondecreasing order of x-coordinates. 

3 Sort P and Q into P y and Q y , respectively, in nondecreasing order of y-coordinates. 

4 S <— ConstructMatrices(Pa;, Q x , P y , Q y , L). 

5 return the rank of t among the elements of sorted matrices in S. 



Figure 3: The ranking algorithm for the Minkowski sum with one linear constraint and a linear 
objective function. 



4 Minkowski Sum Selection with Two Constraints 

In this section, we show the Minkowski Sum Selection problem can be solved in worst-case 
0(n log 2 n) time and expected 0(n log n) time for the case where two linear constraints are given 
and the objective function is linear. 
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4.1 Input Transformation. 



Given P = {(sci,i, yi,i), ■ ■ ■ , (xi >n , yi <n )}, Q = {(^2,1, 2/2,1), (x 2 , n , y2,n)}, a positive integer k, two 
constraints L\\ a\x + b\y > c\ and L 2 : a 2 x + b 2 y > C2, and a linear objective function f(x, y) = dx+ey, 
where 01, b\, c\, a 2 , b 2 , c 2 , d, and e are all constants, we perform the following transformation. 

1. Change the content of P and Q to {(01X1,1 + &ij/i,i, da?i,i + eyi,i), ■ ■ ■ , (ai^i,n + &iyi,n, ^i,n + 
eyi >n )}, and {(aix 2 ,i + 612/2,1, ^2,1 + ey 2 ,i), ■■■ , (aix 2 , n + &i2/2,n, dx 2>n + ey 2 ,n)}, respectively. 

2. Change the constraints from a\x + b\y > c\ and a 2 x + 6 2 y > c 2 to x > ci and + 
Z\l- b b\ a d 2 y ^ c 2> respectively. 

3. Change the objective function from dx + ey to y. 

This transformation can be done in 0(n) time and the answer remains the same. Hence from now 
on, our goal becomes to find the k th largest y-coordinate on the constrained Minkowski sum of P and 
Q subject to the constraints L\\ x > c\ and L 2 : ax + by > c 2 , where a = "llZ^d and b = ^^E^ • 
Note that if the two constraints and the objective function are parallel, we cannot use the above 
transformation. However, if the two constraints are parallel, this problem can be solved in 0(n log n) 
time. For the space limitation, we present the algorithm for this special case in Appendix C. 



4.2 Algorithm 

After applying the above input transformation to our problem instances, there are four possible cases: 
(1) a < 0, b < 0; (2) a > 0, b > 0; (3) a < 0, b > 0; (4) a > 0, b < 0. Note that the two constraints are 
not parallel implies b 7^ 0. If a = 0, we can solve this case more easily in 0(n log 2 n) time by using 
the same technique stated later and we omit the details here. In the following discussion we focus on 
Case (1), and the other three cases can be solved in a similar way. 

For simplicity, we assume that n is a power of two, and each point in {P ®Q)l 1 ,l 2 nas a distinct y- 
coordinate. Now we are ready to describe our algorithm. First, we sort P and Q into P y and Q y , respec- 
tively, in nondecreasing order of y-coordinates using O(nlogn) time. Let P y = {(x^y^), . . . , (x' n ,y' n )} 
and Qy = {(x' lt y[), . . . , (x' n , y' n )}. Denote by Y the sorted matrix of dimensions n x n where the 
(i, j)-th element is y\ + y'j. We then run a loop where an integer interval [I, u] is maintained such that 
the solution is within the set {the h th largest element of Y : I < h < u\. Initially, we set I = 1 and 
u = n 2 . At the beginning of each iteration, we select the "~» +1 -th largest element t of Y, which can 
be done in 0(n) time by Theorem [TJ Let R be the rank of t among the objective values of points in 
(P © Q)l 1 ,l 2 - Then there are three possible cases: (i) R < k; (ii) R = k; (hi) R > k. See Figure H] 
for an illustration. If it is Case (i), then we reset I to u ~ l 2 +1 an d continue the next iteration. If it 
is Case (ii), then we apply the algorithm for the Minkowski Sum Finding problem (discussed in 
Section 6) to find the point p = (x* ,y*) in (P © Q)L!,L 2 ,y<t in 0(n log n) time such that y* is closest 
to t and return y*. If it is Case (iii), then we reset u to and continue the next iteration. 

It remains to describe the subroutine for computing R. Let A = {(x, y)\x < c\ and ax + by < c 2 }, 
B = {(x,y)\x < c\ and ax + by > c 2 and y > t}, C = {(x,y)\x > c\ and ax + by > c 2 and y > t} and 
D = {(x, y)\x > c\ and ax+by < c 2 and y > t}. See Figure[5]for an illustration. First, we compute the 
number of points in (P©Q)n (AuB), say Ri, by calling Rankingi(P, Q, L': x < c\, f'(x, y) = y,t) — 1. 
Secondly, we compute the number of points in (P(BQ) H (AuD), say R 2 , by calling Rankingi(P, Q,L": 
ax + by < c 2 , f'(x, y) = y,t) — 1. Thirdly, we compute the number of points in (P © Q) n A, say P3, 
by calling Rankingi(P, Q, L": ax + by < c 2 , f"(x, y) = —x, c\) — 1. Finally, we compute the number of 
points in (Pffi Q) y >t, say Rt- It can be done by applying Theorem[2]to calculate the rank of t among 
the values of the elements in Y, say R' t , and set R t to R' t — 1. After getting Pi, R 2 , P3, and Rt, we 
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can compute R by the following equation: R = R t — Ri — R 2 + R3 + 1. Since all R%, R2, R3, and R t 
can be computed in O(nlogn) time, the time for computing R is O(nlogra). 




Figure 4: The three possible cases for a given value t. Figure 5: The four regions above 

(exclusive) the line y = t. 

Now let us look at the total time complexity. Since the loop consists of at most O(logn) iterations 
and each iteration takes O(nlogn) time, the total time complexity of the loop is 0(n log 2 n). By 
combining this with the time for the input transformation and sorting, we have the following theorem. 

Theorem 6: Given two n-point multisets PCM 2 and Q C R 2 , a positive integer k, two linear 
constraints L\ and and a linear objective function / : R 2 — » R, the k th largest objective value 
among all objective values of points in (P © Q)l 1: l 2 can be found in 0(n log 2 n) time. 

Theorem 7: For linear objective functions, the Minkowski Sum Selection problem with two 
linear constraints can be solved in expected 0(n log n) time. 

Proof: Due to the space limitation, we leave the proof to Appendix D. □ 



5 Minkowski Sum Selection with A > 2 Constraints 

Let A be a fixed integer greater than two. In the following theorem, we summarize our results of 
the Minkowski Sum Selection for the case where A linear constraints are given and the objective 
function is linear. Due to the space limitation, we leave the proof to Appendix E. 

Theorem 8: Let A be any fixed integer larger than two. The MINKOWSKI Sum Selection problem 
with A constraints and a linear objective function is asymptotically equivalent to the Minkowski Sum 
Selection problem with two linear constraints and a linear objective function. 

6 Minkowski Sum Finding 

In the Minkowski Sum Finding problem, given two n-point multisets P, Q, a set of A inequalities 
Ar > b, an objective function f(x, y) and a real number 5, we are required to find a point v* = (x*,y*) 
among all points in (P © Q)Ar>b which minimizes \ f(x*,y*) — 5\. In this section, we show how to 
cope with an objective function of the form f(x, y) = ax + by or f(x, y) = -| based on the algorithms 
proposed by Bernholt et al. [5j. Instead of finding the point v* = (x*,y*), we would like to focus 
on computing the value of \f(x*,y*) — S\. The point v* = (x*,y*) can be easily constructed from 
the computed information. Before moving on to the algorithm, let us look at the lower bound of the 
problem. 
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Lemma 2: The Minkowski Sum Finding problem with an objective function of the form f(x, y) = 
ax + by or f(x,y) = — has a lower bound of Q(nlogn) in the algebraic decision tree model. 

Proof: Given two real number sets A = {a%, . . . , a n } and B = {bi, ... , b n }, the Set Disjointness 
problem is to determine whether or not A n B = 0. It is known the Set Disjointness prob- 
lem has a lower bound of S7(nlogra) in the algebraic decision tree model [2]. We first prove that 
the Set Disjointness problem is linear-time reducible to the Minkowski Sum Finding problem 
with the objective function f(x, y) = y. Let P = {(1, a%), . . . , (1, a n )}, Q = {(1, — b\), . . . , (1, — b n )}, 
and (x*,y*) be the point in (P © Q) such that y* is closest to zero. Then (x*,y*) = (2,0) if and 
only if A n B ^ 0. Similarly, we can prove that the Set Disjointness problem is linear-time re- 
ducible to the Minkowski Sum Finding problem with the objective function f(x,y) = f. Let 
P = {(1, a%), . . . , (1, a n )}, Q = {(1, —bi), . . . , (1, — b n )}, and (x*, y*) be the point in (P © Q) such that 
£ is closest to zero. Then (x*, y*) = (2, 0) if and only if A n £ ^ 0. □ 

Now, let us look at how to cope with a linear objective function /(x, y) = ax + 6y. Without loss of 
generality, we assume 5 = 0; otherwise we may perform some input transformations first. Thus, the 
goal is to compute the value of min{|ax + by\ : (x, y) G (P © Q)Ar>b}- 

Lemma 3: Divide the xy-plane into two parts: D\ = {(x,y) : ax+by > 0} and D2 = {(x, y) : ax+by < 
0}. Given two points v\ = (xi,y\) and V2 = (X2, y-i) in the same part, let u 7 = (x 7 , y 7 ) = 7^1 + (1— 7)^2, 
where 7 E [0, 1]. Then we have |ax 7 + 6y 7 | > min(|axi + byx\, \ax2 + &2/2 1 ) - 

Proof: It is easy to see the lemma holds if b = 0. Without loss of generality, let x\ > X2 and b 7^ 0. We 
only prove the case where both v\ and v-i are in D\ , and the other case can be proved in a similar way. 
Now consider the following two situations: (1) |axi + 6yi| < |ax2 + &2/2| and (2) |axi+6yi| > \ax2+by^\. 
In the first situation, by \ax% + by%\ < \c1x2 + 67/2! > ax i + ^Ul — 0) an d 0^2 + 62/2 > 0, we can derive that 
^(2/2 — Hi) > o(xi — X2). Let f' = (x', y') satisfy ax\ + 6yi = ax' + by' and x' = x 7 = 7x1 + (1 — 7)2:2. 
It follows that y' = f (1 - 7)(xi - x 2 ) + yi. By y 7 = 7yi + (1 - 7)y 2 = (1 - 7)^2 - yi) + yi, we have 
by~/ > by'. Thus, |ax 7 + 6y 7 | > \ax'+by'\ = \axi + by\\. In the second situation, b{y\ — y2) > a{x2 — x{). 
Let v" = (x",y") satisfy 0x2 + fo/2 = ax" + by" and x" = x 7 = 7x1 + (1 — 7)2:2. It follows that 
V" = §7(^2 - 21) + y 2 . By y 7 = 7^ + (1 - 7)y 2 = 7(2/1 - y 2 ) + y2, we have 6y 7 > 6y". Thus, 
|ax 7 + 6y 7 | > \ax" + by"\ = |ax2 + 6J/2I- Therefore, |ax 7 + 6y 7 | > min (|axi + byi\, \ax^ + by2\) if 
axi + by\ > and 0x2 + fo/2 > 0. □ 

Let D\ = {(x,y) : ax + by > 0} and D2 = {(x, y) : ax + by < 0}. Let i?i be the vertices of the 
convex hull of (P © (5)AT->6,ax+6?/>o and R2 be the vertices of the convex hull of (P © Q) Ar>b, ax+by <o- 
By Theorem [3l we can compute -Ri and R2 in 0(X log A + A • n log n) time. Let soli = min{|ax + by\ : 
(x,y) G (P © Q)Ar>6 H Z?i} and S0I2 = min{|ax + by\ : (x,y) G (P © Q)Ar>6 H ^2}- By LemmaEl we 
have soZi = min{|ax + 6y| : (x, y) G Pi} and S0I2 = min{|ax + by\ : (x, y) G P2j-- Note that both the 
sizes of Pi and P2 are bounded above by 0(A • n). Therefore, sol\ and S0I2 are computable in 0(A • n) 
time by examining all points in Pi and P2. Finally, we have the solution is the minimum of soli an d 
50/2- The total time complexity is 0(Alog A + A • nlogn), and we have the following theorem. 

Theorem 9: Let A be any fixed nonnegative integer. The Minkowski Sum Finding problem with 
A constraints and a linear objective function can be solved in optimal O(nlogn) time. 

Next, we see how to cope with an objective function of the form f(x,y) = -j. Without loss of 
generality, we assume 5 = and a = b = 1; otherwise we may perform some input transformations 
first. Thus, the goal is to compute the value of min{||| : (x,y) G {P®Q)Ar>b\- F° r technical reasons, 
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we define ~ = oo if x = 0. A function / : D — ► (R U oo) denned on a convex subset D of R 2 is 
quasiconcave if whenever vi, t>2 £ -D and 7 G [0, 1] then f(j ■ v\ + (1 — 7) ■ U2) > mrn {/(^2)j /C^)}- 

Lemma 4: Let D 1 = {(x,y) G R 2 : x > 0, y > 0}, D 2 = {(x,y) G R 2 : x < 0, y > 0}, L> 3 = {(x,y) G 
R 2 : x < 0,y < 0}, and D A = {(x,y) G R 2 : x > 0, y < 0}. Define function /j : A M by letting 
fi(x,y) = |f I for each i = 1,2,3,4. Then we have function fa is quasiconcave for each i = 1,2,3,4. 

Proof: We only prove that f\ is quasiconcave. The proofs for f 2 , fs, and can be derived in a 
similar way. Let v\ = (xi,yi) G D%, v 2 = (x 2 ,y 2 ) G D\, and x\ > x 2 . Without loss of generality we 
may assume x\ > and t> 7 G" {^1,^2}- Consider the following two cases. 

Case 1: ^ < ||. Let v' = (x' ,y') be the point which satisfies = f 7 an< ^ x' = x^ = 7x1 + (1— 7)2:2. 
By xi > 0, x 2 > 0, and |i < f|, we have y 2 > ffyi- It follows that y' = (7x1 + (1 - 7)^2)^ = 

jyi + (1 - 7)ffyi < 72/1 + (1 - 7)2/2 = 2/ 7 - By < x' = x 7 and < y' < y 7 , we have = |i > 2L = 
£L = IE, I > min/ISil |2m. 

Case 2: ^ > f|. Let u" = (x",y") be the point which satisfies f| = f" and a;" = x 7 = 
7x1 + (1 — 7)0:2 • By xi > and ^- > we have y% > —y 2 . It follows that y" = (7x1 + (1 — 7)0:2)— = 

X\ 3-2 3-2 3^2 

7f^2/2 + (1 - 7)2/2 < 72/1 + (1 - 7)2/2 = 2/ 7 - By < x" = x 7 and < y" < y 7 , we have 
I — 1 = — >^C = — = 1 — 1> min{ I — I , I — I } • □ 

Let Di = {(x,y) G R 2 : x > 0, y > 0}, D 2 = {(x,y) G R 2 : x < 0, y > 0}, D 3 = {(x,y) G R 2 : 
x < 0,y < 0}, and D4 = {(x,y) G R 2 : x > 0, y < 0}. Let iij be the vertices of the convex hull of 
(P © Q)Ar>b H for i = 1,2, 3, 4. By Theorem [3l each i?j is computable in 0(A log A + A • n log n) 
time. Let soli = mm {|f | : (x,y) G (P © Q)Ar>6 H I^} for each i = 1,2,3,4. By Lemma HJ we have 
soli = min{|f I : (x,y) G Ri} for each i. Note that the size of each Ri is bounded above by 0(A + n). 
Therefore, each soli is computable in 0(A + n) time by examining all points in Ri. Finally, we have the 
solution is min^ =1 soli. The total time complexity is 0(A log A + A • n log n), and we have the following 
theorem. 



Theorem 10: Let A be any fixed nonnegative integer. The Minkowski Sum Finding problem with 
A constraints and an objective function of the form f(x,y) = ^ can be solved in optimal O(nlogn) 
time. 
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Appendix A: Applications to the Length-constrained Sum Selection 
Problem 



Given a sequence S = (s%,S2, ■ ■ ■ ,s n ) of n real numbers, and two positive integers I, u with I < u, 
define the length and sum of a segment <S[i,j] = ...,Sj) to be length(i, j) = j — i + 1 and 
sum(i,j) = Y^h=i s hi respectively. A segment is said to be feasible if and only if its length is in [l,u\. 
The Length-Constrained Sum Selection problem is to find the k th largest sum among all sums 
of feasible segments of S. 

When there are no length constraints, i.e., I = 1 and u = n, the Length-Constrained Sum Se- 
lection problem becomes the Sum Selection problem. Bengtsson and Chen [3] first studied the Sum 
Selection problem and gave an 0(n log 2 n)-time algorithm for it. Recently, Lin and Lee provided an 
0{n log n)-time algorithm [13J for the Sum Selection problem and an expected 0{n log(u — I + 1))- 
time randomized algorithm [H] for the Length-Constrained Sum Selection problem. In the 
following, we show how to solve the Length-Constrained Sum Selection problem in worst-case 
0(nlog(u — I + 1)) time. 

Algorithm 

We first reduce the Length-Constrained Sum Selection problem to the Minkowski Sum Selec- 
tion problem as follows. Let P = {po,Pi, ■ ■ ■ ,Pn} and Q = {qo, g_i, . . . , q- n }, where pi = (xi, yi) = 

(*i Y!t=i s t) and Q-i = &i,Vi) = (-«> - Et=i s t) fo r all i = 0, 1, ... ,n. 

A point (x, y) in P(B Q is said to be a feasible point if and only if I < x < u. Each feasible segment 
S[i,j] corresponds to a feasible point (x,y) = pj + qi-i in P ®Q. Thus, the Length-Constrained 
Sum Selection problem is equivalent to finding the k th largest y-coordinate among all y-coordinates 
of feasible points in P © Q. We next show how to do this in 0(n log(n — I + 1)) time. For simplicity, 
we assume n is a multiple of u — I. 

1. Let i t = t(u - I) and j t = I — t(u - I) for t = 0, 1, . . . , 

2. For t <- to ^_ do 

u — t 

(a) Let P t ={ph e P : it - (u-l) < h <i t } and Q t = Q t ,i U Q t ,2, where Q t) i = {qn G Q ■ jt < 
h<jt + (u - I)} and Q t , 2 = Uh G Q : jt + (u - I) < h < j t + 2{u - I)}. 

(b) Store the y-coordinates of points in (P t © Qt,i)x>i as a set N tt \ of sorted matrices such 
that the sum of side lengths of the sorted matrices in iV^i is no greater than c • {{\Pt\ + 
|Q t) i|)log(|Pt| + |<5i,i| + 1)) for some constant c. 

(c) Store the y-coordinates of points in (P t © Qt,2)x<u as a set N t: 2 of sorted matrices such 
that the sum of side lengths of the sorted matrices in Nt2 is no greater than c • {{\Pt\ + 
\Qt,2\) log(|-Pf| + \Qt,2\ + 1)) f° r some constant c. 



3. Return the k th largest element among the elements of sorted matrices in (Jt=o C^t,i ^ ^,2)- 
The following lemma ensures the correctness. 



Lemma 5: (P © Q)i< x < u = U t to° ((-P* © Qt,i)/<* U (P t © Qt,i)x<n). 
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Proof: We prove that (P © Q)i< x < u = U=o° ( p t © Q)l<x<u = Ut=o° © Qt)i<x<u = Ut=o° ((^ © 

Qt,i)KrE<« u { p t © Qt,2)i<x<u) = Ut=o° ((-^ © Qt,i)i<x U (P t © Qt,2)x<«)- It is clear that equations (1) 
and (3) are true, so only equations (2) and (4) remain to be proved. 

We first prove equation (2) by showing that (Pt © Q)i< x <u = (Pt © Qt)i<x<u- Suppose for contra- 
diction that there exist pi £ P t and qj Q t such that I < X{ + Xj = i + j < u. By pi G Pt, we have 
(t — l)(u — I) < i < t(u — I); by qj Qt, we have either j < I — t(u — I) or j > I — (t — 2)(u — I). 
It follows that i + j is either less than I or larger than u, a contradiction. To prove equation (4), 
it suffices to prove that all points in (Pt © Qti) must have ^-coordinates less than u and all points 
in (P t © Qt,2) must have x-coordinates larger than I. Let pi £ P t , qj £ Qt,i and qj> S Qt,2- It 
follows that (t — l)(it — I) < Xi = i < t(u — I), I — t(u — I) < x j = j < I — (t — l)(u — I), and 
l — (t — \)(u—l) < Xj' = j' < l—(t — 2)(u—l). Thus, we have x-i+Xj = i+j < u and / < x-i+Xji = i+f. □ 

Since \Pt\, |Qt,i|, and \Qt,2\ are no greater than (u — I) for all t, each execution of Step 2.b and 
Step 2.c can be done in 0((u — I) log(u — I + 1)) time by LemmaHJ There are total ^37 + 1 iterations of 
the for-loop in Step 2, so the total time spent on Step 2.b and Step 2.c is 0(n\og(u — l + l)). The sum of 

n, n 

side lengths of sorted matrices in (J^ (N t ,iUN t , 2 ) is 0(Z t =o ELiG^I + IQmI) log(|P| + |Q M | + l)) = 

n 

O(j2t=o( u ~ — ' + 1)) = 0(n\og(u — I + 1)). Therefore, by Theorem[[J Step 3 can be done 

in 0(nlog(u — 1 + 1)) time. Putting everything together, we have that the total running time is 
0(n\og(u - I + 1)). 

Theorem 11: The Length-Constrained Sum Selection problem can be solved in 0(n \og(u — 
I + 1)) time. 

Appendix B: Applications to the Density Finding Problem 

Given a sequence of number pairs S = ((s\,wi), (^2,^2), • • • , (sm %)) where Wi > for i = 1, 2, . . . ,n, 
two positive numbers l,u with / < u, and a real number S, let segment S(i,j) of S be the consecutive 
subsequence of S between indices i and j. Define the sum s(i,j), width w(i,j), and density d(i,j) 
of segment S(i,j) to be J2i=i s r, J2t=i w r and ^ferjj ; respectively. A segment S(i,j) is said to be 
feasible if and only if / < w(i,j) < u. The Density Finding problem is to compute the density of 
the feasible segment S(i* ,j*) which minimizes \d(i*,j*) — S\. Lee et al. [12] proved that the Density 
Finding problem has a lower bound of £l(n\ogn) in the algebraic decision tree model and provided 
an 0(n log 2 m) algorithm for it, where m = min( I u ~ l I , n) and w m i n = mini< r < n w T . In the following 
we describe how to solve the Density Finding problem in 0(n log n) time by using the algorithm 
developed in Section El 

Let io(l,0) = and s(l,0) = 0. Compute in 0(n) time the following two point sets: P = 
{(w(l,i), s(l,i))\0 < i < n} and Q = {(— w(l, i), — s(l, i))\0 < i < n)}. Note that each feasible 
segment S(i, j) of S corresponds to a point (w(l, j) — w(l, i — 1), s(l,j) — s(l, i — 1)) in (P © Q)i< x < u - 
Thus, the problem is reduced to finding the point (x*,y*) in (P © Q)i< x <u such that \^ — S\ is 
minimized. By Theorem 1101 it can be done in 0(n log n) time, so we have the following theorem. 

Theorem 12: The Density Finding problem can be solved in optimal 0(n log n) time. 

Appendix C: Minkowski Sum Selection with Two Parallel Constraints 

Now we explain how to solve the Minkowski Sum Selection problem with two parallel constraints in 
O(nlogn) time. Given P = {(x ltl ,y lt i), (x 1>n ,y 1>n )}, Q = {(x 2) i, 2/2,1), • • • , (x2,n,y2,n)}, a positive 



14 



integer k, two parallel constraints L%: ax + by > c\ and L 2 : ax + by < c 2 with c\ < c 2 , and a linear 
objective function f(x,y) = dx + ey, where a, 6, ci, C2, d, and e are all constants, we want to find 
the k th largest objective value among all objective values of points in (P © Q)li,l 2 - Note that if the 
constraints L\ and L 2 are of the forms ax + by < c\ and ax + by < C2 respectively, this problem 
degenerates to the Minkowski Sum Selection problem with one constraint and can be solved by 
the algorithm stated in Section 3. We perform the following transformation. 

1. Change the content of P and Q to {(0x1,1 + 6y lj i, 0x1,1 + eyi ) i), . . . , {ax\^ n + 6yi, n , dxi tTl + eyi jn )}, 
and {(0x2,1 + 6y 2 ,i, dx 2 ,i + ey 2 ,i), • • • , (ax 2 , n + by 2 , n , dx 2 , n + ey 2) n)}, respectively. 

2. Change the constraints from ax + by > c\ and ax + by < c 2 to x > c\ and x < c 2 , respectively. 

3. Change the objective function from dx + ey to y. 

This transformation can be done in 0(n) time and the answer remains the same. Hence from now 
on, our goal becomes to find the k th largest y-coordinate on the constrained Minkowski sum of P and 
Q subject to the constraints L\: x > c\ and L 2 : x < c 2 . First we sort P and Q into P x and Q x in 
nondecreasing order of ^-coordinates, respectively in O(nlogn) time. Let P x = {(xi,y%), . . . , (x n , y n )} 
and Q x = {(xi,y~i), • • • , (x n ,y n )}. For all points in P x , we can form a partition of them according to 
the values of their ^-coordinates. Let the partition be Pt i; Pt 2 , . . . ,Pt m where Pt t = {(xj,yj) £ P x : 
{fi — l)(c2 — ci) < Xj < U(c 2 — ci)}, and t{ be an integer for i = 1, 2, . . . , m with t\ < t 2 < ■ ■ ■ < t m . 
Similarly, we can partition the points in Q x according to the values of their x-coordinates. Let the 
partition be Q fl , Q t - 2 Qf m , where Q f . = {(xj,yj) G Q x : ci-ij(c 2 -ci) < xj < ci-(ti-l)(c 2 -ci)}, 
and ti be an integer for i = 1,2, ... ,m' with t\ < t 2 < ... < t m >. Since P x and Q x are sorted in 
nondecreasing order of x-coordinates respectively, the two partitions can be easily produced in linear 
time. In the following, we show the algorithm for this problem. 

1. Let Pt x , Pt 2 , • • • , Pt m and , Qi 2 , . . . , Qi m , be defined as the above. 

2. For i <— 1 to m do 

(a) To find Q u = {(xj,y~j) £ Q x : c\ - U(c 2 - c\) < Xj < c\ - (U - l)(c 2 - ci)} and Q tl -\ = 
{{xj,y~j) G Q x : ci - (U - l)(c 2 - ci) < % < ci - (U - 2)(c 2 - ci)}. 

(b) If Qt 4 exists, store the y-coordinates of points in (P t . © Qtjx^ci as a set iV^i of sorted 
matrices such that the sum of side lengths of the sorted matrices in JVj 1 is no greater than 
c ' (0-ftJ + logdi^J + |QtJ + 1)) for some constant c. 

(c) If Qti-i exists, store the y-coordinates of points in (P ti © Qti-i)x<c 2 as a set Ni j2 of sorted 
matrices such that the sum of side lengths of the sorted matrices in Ni t2 is no greater than 
c ' ((\Pti\ + \Qu-i\)\og(\P tl \ + \Qti-l\ + 1)) for some constant c. 

3. Return the k th largest element among the elements of sorted matrices in U^LiO^j.l ^ -^,2)- 

By the proof of Lemma \5\ we ensure the correctness of the algorithm. Now we consider the time 
complexity of the algorithm. By Lemma [H each execution of Step 2.b and Step 2.c can be done in 
C((|P ti | + |Q ti |)log(|P ti | + |g ti | + l)) andOaiPj + lQ^DlogdPJ + lQ^I + l)) time, respectively. 
Since there are total m iterations of the for- loop in Step 2, it follows that 

m 

E (d p **l + I^DlogdPj + \Q tl \ + 1) + (IPJ + |Q ti _i|)log(|P ti | + |Q ti _!| + 1)) 
i=i 
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m 

< E (0 P **I + \Qu\)^g(2n + 1) + (\P U \ + IQ^I) log(2n + 1)) 

m 

< ^(2|P t .| + |Q ti | + |Q ti _ 1 |)log(2n + l) < (4n) log(2n + 1). 
i=l 

Therefore, the total time spent on Step 2 is O(nlogn) and the sum of the side lengths of sorted 
matrices in Ui=i(-^j>i UJV^) is also O(nlogn). By Theorem[IJ Step 3 can be done in O(nlogn) time. 
Putting everything together, we have that the total running time is 0(n log n). The next theorem 
summarizes the time complexity of the algorithm. 

Theorem 13: For linear objective functions, the Minkowski Sum Selection problem with two 
parallel constraints can be solved in O(nlogn) time. 

Appendix D: A Randomized Algorithm for Minkowski Sum Selection 
with Two Constraints 

In this section, we introduce a randomized algorithm for the Minkowski Sum Selection problem 
with two constraints that runs in expected 0(n log n) time. 

Subroutines for Minkowski Sum Selection Problem with Two Constraints 

Our randomized algorithm is based on three subroutines for three subproblems. In this subsection, 
we define these subproblems and give these subroutines for them. 

Before we discuss these subroutines, we introduce the notion of an order-statistic tree. An order- 
statistic tree is a balanced search tree with additional information, size[z], stored in each node z of the 
tree. The additional information size[z] contains the total number of nodes in the subtree rooted at z. 
Define left[z] and right[z] are the left and right children of the node z, respectively. The additional 
information size[z] equals to size[left[z]] + size[right[z]] + 1 if z is an internal node, and one if z is 
a leaf node. Let key[z] be the key of the node z. The rank of a given value x can be determined in 
O(logn) time by using the order-statistic tree T, where n is the number of nodes in T. That is, we 
can find the rank r(x,T) = \{y\y £ T,key[y] > x}\ in O(logn) time, retrieve an element in T with 
a given rank in O(logn) time and maintain both insertion and deletion operations in T in O(logn) 
time. 

The first subproblem is the reporting version of the Minkowski Sum Range Query problem 
with two constraints, which is defined as follows: Given two n-point multisets P, Q, two constraints 
Li, L2, and two real numbers s/, s r with s; < s r , we want to output all points in (P © Q)l 1 ,l 2 i such 
that their ^-coordinates are in the range [si, s r ]. Before we discuss this subproblem, we consider a weak 
version of this subproblem. The weak version is defined as above, except that we assign the constraints 
L\\ ax + by > c\ and Li'- ax + by < C2, where c\ < oi. To solve the weak version, we perform the 
input transformation stated in Section 4.1.1 to change L\ and L2 to x > c\ and x < C2 respectively. 
Then we sort P and Q into P x and Q x respectively in nondecreasing order of x-coordinates. Let 
P x = ((xi,yi), • • • , (x n ,y n )) and Q x = ({xi,yi), ... , (x n ,y n )). Denote by k the smallest index in Q x 
such that Xi + > c\ and by r , the largest index in Q x such that x% + x Ti < 02- It is guaranteed 
that li > lj and rj > rj for any i < j. For this reason, it can be easily done to find all li and ri for 
i = 1, . . . , n in total 0(n) time. We use the points in Q x to construct an order-statistic tree with the 
values of y-coordinates as keys. For the index j, we use the points (xi,yi) in Q x for i = lj,. . . ,rj to 
construct an order-statistic tree T(j) in 0((rj — lj) log(rj — lj)) time. Because T(j) is also a balanced 
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binary search tree, we can report all points in T(j) whose y-coordinates are in [s; — yj,s r — yj] by 
binary search in 0(log(rj — lj) + hj) time, where hj is the total number of points whose y-coordinates 
are in [sj — yj, s r — yj] . At each iteration j, we can maintain T(j + 1) dynamically by deleting all points 
(xi, y~i) in T(j) for i = rj + i + 1, . . . , Tj and inserting all points {xy , y>) into T{j) for i' = ... — 
It suffices to iterate on each index j to find all points in (P © Q)L 1 ,L 2 ,y<s r ,y>s r Hence, we can solve 
the weak version in 0(n log n + h) time, where h is the output size. 

Now we show how to solve the reporting version of the Minkowski Sum Range Query problem 
with two constraints by the weak version stated above. By performing the transformation stated in 
Section 4.1.1, L% and L<i are changed to x > c\ and ax + by > C2 respectively. We divide the feasible 
region bounded by L\, L-2, and s; < y < s r into several subregions. For each subregion, it can be 
solved by the subroutine for the weak version of this subproblem. Let L\ be the line that is parallel 
to L\ and passes through the intersection point of L2 and y = s r , and L 2 be the line that is parallel 
to L2 and passes through the intersection point of L\ and y = s r . By the location of the intersection 
point of L\ and L' 2 , we have four possible cases: (a) the intersection point lies in the line y = s\; (b) 
the intersection point lies below the line y = si\ (c) the intersection point lies above the line y = s r ; 
(d) the intersection point lies between y = S[ and y = s r . See Figure [6] for an illustration. We solve 
the reporting version according to the four possible cases respectively. For Case (a), we consider the 
parallelogram formed by L2, L' 2 , y < s r , and y > si and all feasible points in this parallelogram can 
be reported by the subroutine of the weak version. Next we consider the rectangle formed by L\, L\, 
y < s r , and y > si and all feasible points in this rectangle can be reported in the same way, except we 
have to remove the redundant points in the area formed by Uy, L' 2 , and y < s r . When we report each 
feasible point in the rectangle formed by L±, L\, y < s r , and y > si, we can check whether this point 
lies in the area formed by L^, L 2 , and y < s r in O(l) time. If this point lies in this area, we discard 
this point, or we report it. For Case (b), we can also solve this case in the same way of Case (a), 
except the redundant points are in the area formed by L' l5 L' 2 , y < s r , and y > si. For each redundant 
point, however, the removal can be also easily done in 0(1) time. 

For Case (c), we must divide the triangle formed by L\, L2, and y > si in another way. Let a be 
the intersection point of L\ and L2, c be the intersection point of L\ and y > si, and b be the middle 
point of the line segment oc. We can draw the line L3 that is parallel to y = si and passes through b, 
and let d be the intersection point of L3 and L2. Let L 2 be the line that is parallel to L2 and passes 
through b, and L" be the line that is parallel to L\ and passes through d. Because the triangle formed 
by Li, L2, and y > si is a right-angled triangle, the intersection point of L'{ and L 2 must lie in the 
line y = si. We first report the feasible points in the parallelogram formed by L%, L", L2, and L 2 by 
the subroutine of the weak version. Then we report the feasible points in the parallelogram formed 
by L3, y > si, L2, and L 2 in the same way and remove the redundant points, i.e., the points lie in the 
area formed by L3, L", and L 2 . Finally, we report the feasible points in the rectangle formed by L\, 
L'(, L3, and y > si and remove the redundant points in the area formed by L3, and L 2 . 

For Case (d), we report all feasible points in the parallelogram formed by L2, L 2 , y < s r , and 
y > S[ and in the rectangle formed by L\, L^, y < s r , and y > si with the removal of the redundant 
points in the area formed by L\, L 2 , and y < s r . The remaining is the triangle formed by L\, L' 2 , and 
y > S[ and is just Case (c). Therefore, we can solve this triangle in the same way of Case (c). In the 
following, we conclude the time complexity of this problem. 

Lemma 6: The reporting version of the Minkowski Sum Range Query problem with two con- 
straints can be solved in 0(n log n + h) time, where h is the output size. 

The second subproblem, called the counting version of the Minkowski Sum Range Query 
problem with two constraints, is defined as before, but we only want to find the number of points 
in (P @ Q)l!,l 2 satisfying the range query, i.e., their y-coordinates are between s; and s r . To solve 
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Figure 6: The illustration of how to solve the reporting version. 



this subproblem, we make use of the procedure Ranking2 shown in Figure [7J Let Ri be the number 
of points in (P © Q)l ± ,l 2 that their ^-coordinates are larger than si. It is obvious to see that Ri 
is equal to Ranking2(P, Q, L%, L2, f(x, y) = y,si). Let R r be the number of points in (P © Q)l 1 ,l 2 
that their y-coordinates are larger than s r , and R r is equal to Ranking2(P, Q, L\, L2, f(x, y) = y, s r ). 
As a result, the number of points in (P © Q)L 1 ,L 2 ,y<s r ,y>s l is Ri — R r phis the number of points in 
(P®Q)L 1 ,L 2 ,y=s r Let R be the number of points in {P(BQ)L 1 ,L 2 ,y=s r R can be obtained by computing 
Rankingi(P, Q,y = s/ ; f(x, y) = x, c\) minus Rankingi(P, Q,y = si, f(x, y) = x, C2). Hence, we obtain 
the number of points in (P © Q)L 1 ,L 2 ,y<s r , y >s r 

Lemma 7: The counting version of the Minkowski Sum Range Query problem with two con- 
straints can be solved in 0(n log n) time. 

The last subproblem, called the Random Sampling Minkowski Sum problem with two con- 
straints, is defined as follows: Give two n-point multisets P, Q, two constraints L\\ x > ci, L2: 
ax + by > C2, and two real numbers Si, s r with si < s r , we want to randomly generate n points from 
L 2 ,y<s r ,y>si with replacement. For the ease of similar discussions, in the following we only 
focus on Case (a) illustrated in Figure [6] and the other cases can be solved in a similar way. 

Let iV be the number of points in (P © Q)L 1 ,L 2 ,y<s r ,y>s v ^i be the number of points in the 
parallelogram A formed by L2, L' 2 , y < s r and y > si, N2 be the number of points in the rectangle B 
formed by L\, L' x , y < s r and y > si, and N3 be the number of points in the triangle C formed by L\ , 
L' 2 and y < s r . N, Nx, N2, and A^3 can be computed by the subroutine for the counting version of 
the Minkowski Sum Range Query problem with two constraints. For the parallelogram A, we can 
use the subroutine for the weak version of the reporting subproblem to construct the order-statistic 
tree T(j) on points in Q x for j = 1, . . . , n. Let the size of T(j) be tj and then N\ = t\ + 12 + ■ ■ ■ + t n . 
We first pick n random integers R = {f\,f2, ■ ■ ■ ,r n } uniformly distributed in the range from 1 to 
./V with replacement. Since N = 0(n 2 ), we can sort them by radix sort and rename them such that 
r\ < T2 < ■ ■ ■ < f n m 0(n) time. Let Tj = t\ -K2 + • • •+£/■ For each j, there exist f c , f c+ \, . . . , f c+ d 6 R 
such that Tj-i < f c < f c+ \ < . . . < f c+c i < Tj. For fj < JVi, we shall find points s c , s c+ ±, . . . , s c+ d 
in (P © Q)Li,L 2 ,y<s r ,y>si with a one-to-one correspondence to f c ,f c+ i, . . . ,f c+( [. For each index j, we 
make a query to the order-statistic tree T(j) in order to count the total number aj of points such 
that their y-coordinates are less than s/, i.e., aj = \{y\y G T(j),key[y] < si}\. We then retrieve the 
point s c+ i = (x qi ,y qi ) from T(j) such that y qi has a rank equal to aj + f c+ j — Tj-_i in T(J) for each 
i = 0, 1, . . . , d. For fj > N\, we record the total number of fj larger than N\, say (j). We thus obtain a 
set of points, S = {s\, . . . , s^}. Next we remove the points in S lying in the triangle C, and record 
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the number of points removed, say tp. Then we select (f) + tp points randomly from the rectangle B 
using the same way. Combining these 4> + tp points with S, we obtain the set of random sampling 
points. 

Now we show that the sampling resulted from our random sampling subroutine is uniformly ran- 
dom. Let X = {Xi,X 2 , ■ ■ ■ ,X n } be any fixed random sample generated from our random sampling 
subroutine and X^ , X^ 2 , . . . Xi k be the points in the triangle formed by L\, L 2 , and y > sj. The 
probability that this random sample occurs is exactly (j[) k (jf) n ~ k (^) n ^ k = (jf) n - Thus we obtain 
that for any fixed random sample, the probability of occurring is the same (jf) n , he., the random 
sample generated from our random sampling subroutine is uniformly random. The following lemma 
concludes the time complexity. 

Lemma 8: The Random Sampling Minkowski Sum problem with two constraints can be solved 
in expected (9(n log n) time. 

Algorithm 

First of all, we perform the input transformation stated in Section 4.1.1. Hence we assume that the 
two constraints are L\\ x > c±, L 2 : ax + by > c 2 and the objective function is f(x,y) = y. We say a 
range [si, s r ] contains a point s £ P © Q if s satisfies two linear constraints y>s\ and y < s r . 

Let the y-coordinates of the points in P®Q be in the range [si, s r ], and s* = (x*,y*) be the solution 
of the Minkowski Sum Selection problem with two constraints. Our randomized algorithm for this 
problem will contract the range [s;,s r ] into a smaller range [s^,s r '] such that [s;/,s r '] contains s* . A 
point s is said to be feasible if s € (P© Q)l 1 ,l 2 - Let N be the number of points in (P © Q)l 1 ,l 2 - The 
subrange [sf, v] will contain at most 0(N/ y/n) feasible points. We shall repeat the contraction several 
times until the subrange contains at most 0(n) feasible points and also the solution s* . Then we output 
all feasible points in this subrange by the subroutine for the reporting version of the Minkowski Sum 
Range Query problem with two constraints and find the solution s* with an appropriate rank. 

We first randomly select n independent feasible points S = {si, s 2 , ■ ■ ■ , s n } which are contained in 
the range [si,s r ] from TV feasible points by the subroutine for the Random Sampling Minkowski 
Sum problem in O(nlogn) time. When we randomly select a feasible point in [si, s r ], the probability 
is jj that its y-coordinate is smaller than that of s*. Consider this event as a Bernoulli trial with 
the success probability p = jj. It is obvious to see that the total number of successes is a random 
variable which has a binomial distribution. Hence the expected value of the total number of successes 
is fi = np = n^f. As a result, we know the good approximation for the point of the k th largest 
y-coordinate among all feasible points is the point of the e th largest y-coordinate among S, where 
e = [np\ = [n^\ . 

Let I' = max{l, [n-j^ — t^-\} and r' = min{n, [n^ + i^J}, where t is a constant and will be 
determined later. After the random sampling, we can find the I'-th and the r'-th largest y-coordinates 
in S, say s// and s r i respectively, by any standard selection algorithm in 0{n) time to obtain the 
subrange [si>,s r i]. Next, we check the following two conditions by the subroutine for the counting 
version of the Minkowski Sum Range Query problem with two constraints in O(nlogn) time: 

(1) The subrange [sj/, v] contains the solution s*. 

(2) The subrange [sj/,s r /] contains at most t 2 N/(t — l)^/n (< 2tN/yJn) feasible points. 

Let k\ and k 2 be the total number of feasible points contained in [si, s^) and [si, s r /] respectively. If 
s* is contained in the subrange [s^, s r >], we know that k± < k and k 2 > k. If both of the two conditions 
hold, we replace the range [si, s r ] and the rank k with the subrange [s^, s r i] and the rank k' = k — k±. 
If any of the two conditions is violated, we repeat the above step until both of the two conditions 
are satisfied, i.e. we need to select n random feasible points with replacement in the range [s;,s r ] by 
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running the subroutine for the Random Sampling Minkowski Sum problem again to obtain a new 
subrange [s[i, s r i] and then check the above two conditions for this new subrange. 

Since this randomized algorithm for the Minkowski Sum Selection problem with two con- 
straints starts with TV" feasible points, after the first successful contraction, we have a new range 
[s;/, s r r] contains 0(N/ \/n) feasible points and s* , the point of the k'-th largest y-coordinate among all 
feasible points. After the second successful contraction, we have an another subrange [sj«,s r «] which 
contains 0( N ^ ) = 0(n) feasible points and s* , the point of the k"-th largest y-coordinate among 
all feasible points. Since the number of feasible points contained in the range [si",s r n] is 0(n), we 
can enumerate all feasible points in this range in 0(n log n) time via the subroutine for the reporting 
version of the Minkowski Sum Range Query problem with two constraints and select the point of 
the k"-th largest y-coordinate from these feasible points by using any standard selection algorithm in 
0(n) time. 

Now we show that with high probability, the point of the k th largest y-coordinate among all feasible 
points contained in subrange [sir, s r i] and the subrange [s//, s r >] contains at most t 2 N/(t—l)^/n feasible 
points. Applying the results of Matousek et al. [19], we have the following lemma: 

Lemma 9: Given a set of feasible points = {61,82, ■ ■ ■ ,8n}, an index k (1 < k < N), and an 
integer n > 0, we can compute in O(nlogn) time an interval [s^,s r /], such that, with probability 
1 — l/Q,{y/n), the point of the k th largest y-coordinate among lies within this interval, and the 
number of points in that lie within the interval is at most N/£l(y/n). 

By the results given by Matousek et al. |19j . we can choose t = 3. Therefore, we just need to 
repeat the contraction step at most twice on average in the randomized algorithm for the Minkowski 
Sum Selection problem with two constraints. We conclude the time complexity of the randomized 
algorithm in the following theorem. 

Theorem 14: For linear objective functions, the Minkowski Sum Selection problem with two 
linear constraints can be solved in expected 0(n log n) time. 

7 Appendix E: Minkowski Sum Selection with A > 2 Constraints 

In the following, we describe how to solve the Minkowski Sum Selection problem with A > 2 
constraints and a linear objective function by using Algorithm Selection and Algorithm Ranking2. 
The pseudocodes of Algorithm Selection and Algorithm Ranking2 are given in Figure [7] and Figure [HJ 
respectively. The running time of Algorithm Selection is 0{n log 2 n) by Theorem El and the running 
time of Algorithm Ranking2 is 0(n log n). 

Without loss of generality, we assume the given objective function is f(x, y) = y; otherwise, we 
may perform some transformations on the input first. Before moving on to the algorithm, let us 
pause here to introduce some definitions. Let x = {Li,L,2, . . . ,L\} the set of the A constraints and 
V = (v\ = (x\,yi), V2 = (x2,y2), ■ ■ ■ , v m = (x m , y m )) be the vertices of the polygon formed by the A 
constraints, sorted in nonincreasing order of y-coordinates. For simplicity, we assume that each point 
in (P © Q)x an d V has a distinct y-coordinate and the polygon formed by the A constraints is closed. 
If not, we can add another constraint making that the polygon is closed and contains all points in 
(P © Q)x- Denote by L\ and L l 2 the two constraints resulting in the edges of the polygon between 
lines y = yi and y = yi + \ for each i = 1, . . . , m — 1. For each i = 1, 2, . . . , m — 1, we define Ri to be 
the point set {(x,y) : (x,y) G (Pffi Q) x and yj + i < y < yi}. An illustration of the above definitions is 
shown in Figure [H 
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Algorithm Selection (P, Q, L\, L 2 , f, k) 

Input: PCM 2 and QCR 2 are two multisets; L\ and L2 are linear 

constraints; / : IR 2 — ► R is a linear objective function; k is a positive integer. 
Output: The k th largest value among all objective values of points in (P © Q)li,l 2 - 

1 Perform the input transformation in Section 4.1. Let L\ : x > C\ and 

L 2 : ax + by > C2 be the resulting constraints after the input transformation. 
/*Assume that a < and b < 0.*/ 

2 Sort P and Q into P y and Q y , respectively, in nondecreasing order of y-coordinates. 

3 Let P y = ((xi,2/i),...,(a4,?4)) and Q y = ((x[, y[), . . . , (x' n , y' n )). 

4 Denote by Y the sorted matrix of dimensions n x n where the (i,j')-th element 
is y'i + y'j- 

5 u <— n 2 ; I <— 1; m <— ; t <— the m i/l largest element of 1". 

6 while true do 

7 else 



8 i?<- Ranking 2 (P,Q,Li,L 2 ,/,t). 

9 if R < k then 

10 I <— m; m <— " + 2 +1 ; t <— the m" 1 largest element of Y. 

11 else if R = k then 

12 Find the point p = (x*,y*) in (P © Q)L li L 2 ,x<t such that y* 

is closest to t. 

13 return y*. 

14 else 

15 u <— m; m <— ; t <— the largest element of Y. 



Figure 7: The algorithm for the Minkowski Sum Selection problem with two linear constraints 
and a linear objective function. 

We now explain how our algorithm works. First of all, we have to compute the value of |Pj| for 
each i = 1, 2, . . . , m — 1. The value of \Ri\ is equal to the number of points in (P © Q)^ £» above the 
line y = y^+i minus the number of points in (P © Q)^^ above the line y = y^. The number of points 
in {_P®Q)l\,l\ above the line y = yj+i is equal to the rank of yj+i among the values of y-coordinates of 
the points in (P®Q) L ^ L » minus one, i.e., Ranking2(P, Q, L\,L % 2 , f(x, y) = y, yj+i) — 1. The number of 
points in (P©Q) L i L » above the line y = yi is equal to the rank of yi among the values of y-coordinates 
of the points in (P © Q)l\,l\ mm us one, i.e., Ranking 2 (P, Q, L\, L 2 , f(x, y) = y,y«) — 1. Thus, we 
have \Ri\ = Ranking 2 (P, Q, L\,L l 2 , f(x, y) = y, y i+1 ) - Ranking 2 (P, Q, L\,L\, f(x, y) = y, yi) for each 
i = 1, 2, . . . , m — 1. After computing all values of |Pj| for each i = 1, 2, . . . , m — 1, we can determine 
in which set P w the solution y* is located. 

Next, we explain how to determine this set R w where the solution y* is located. Let Sq = and 
Si = \R\ I + |P 2 | + • • • + \Ri\ for each i = 1, 2, . . . , m— 1. We find the smallest index u; such that S w > k 
or w = m — 1. It is easy to see that the solution y* must be located in the set Rw. It follows that the 
solution y* must be among the y-coordinates of points in (P © Q)l™,l%- 

The last step is to obtain the rank of y* among the y-coordinates of points in (P © Q)l™,l™ 
and then we can use Algorithm Selection 2 to find the solution y*. Because the solution y* is lo- 
cated in the set R w , the rank of y* among the y-coordinates of points in R w is k — S w -i. We call 
Ranking 2 (P, Q, Lf, L™, f(x, y) = y,y w ) and let r be the return value of the invocation. The value 
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Algorithm Ranking2 (P, Q, L\, L2, /, t) 

Input: Two multisets PCM 2 and Q C R 2 ; two linear constraints L\ and 

a linear objective function / : R 2 — > R; a real number t. 
Output: The rank of t among the objective values of points in (P © Q)li,l 2 - 

1 Perform the input transformation in Section 4.1. Let Li : a; > ci and 

L2 : ax + by > c 2 be the resulting constraints after the input transformation. 
/* Assume that a < and b < 0.*/ 

2 Sort P and Q into P y and Q^, respectively, in nondecreasing order of y-coordinates. 

3 Let P y = {(x' l ,y' l ),...,(x' n ,y' n )) and Q y = ({x[, y[), . . . , (x' n , y' n )). 

4 Denote by Y the sorted matrix of dimensions n x n where the )-th element 
is Vi + y'j- 

5 fii <- Rankingi(P, Q, L': x < c\,f'(x,y) = y,t) — 1. 

6 R2 <— Rankingi(P, Q, L": ax + by < C2, f'(x, y) = y,t) — 1. 

7 P 3 <— Rankingi(P, Q, L": ax + by < C2,f"(x,y) = -x,ci) - 1. 

8 Pj <— the rank of i among the values of y-coordinates of the points in K minus one. 

9 R <- R t - Pi - P 2 + P 3 + 1. 

10 return P. 



Figure 8: The ranking algorithm for the Minkowski sum with two linear constraints and a linear 
objective function. 

r — 1 is equal to the number of points in (P © Q)l^,l% that are above the line y = y w . Hence, 
there are r points in (P © Q)l™,l™ that are above the line y = y w and there are k — S w -i points 
in (P © Q)l™,l™ belonging to the set R w . It is easy to derive that the value r — 1 + k — S w -± is 
the rank of the solution y* among the y-coordinates of points in (P © Q)l™,l^- Finally, we call 
Selection2(P, Q, Lf, L™, f(x, y) = y,r — 1 + k — S w -i) to find the solution y*. The detailed algorithm 
is given in Figure [TUl 

Now let us see the running time. Because there are A constraints, it takes O(AlogA) time to 
compute V [4]. To compute the values of |Pj| for all % G {1,2,... , m — 1}, we call Algorithm Ranking2 
2A times, which takes 0(A • ralogn) time. The invocation of Algorithm Selection takes 0(nlog 2 n) 
time by Theorem [H Hence, the total time complexity is 0(X ■ nlogn + A log A + nlog 2 n). Note that 
if we use the randomized algorithm described in Theorem [7] instead of Algorithm Selection, the total 
time complexity is 0(A • n log n + A log A + n log n). Hence, we have the following theorem. 

Theorem 15: Let A be any fixed integer larger than two. The Minkowski Sum Selection problem 
with A constraints and a linear objective function is asymptotically equivalent to the Minkowski Sum 
Selection problem with two linear constraints and a linear objective function. 
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Figure 9: The polygon formed by the A constraints. 



Algorithm Selection a (P, Q, x, k) 

Input: Two multisets P C IR 2 , Q C M 2 ; a set x of the A constraints; a positive integer 
Output: The A; i/l largest value of y-coordinates among all points in (P © Q)l 1 ,...,L\- 

1 V <— vertices of the polygon formed by the A constraints in x- 

2 m <- |V|; So <- 0. 

3 for i <— 1 to m — 1 do 

4 <- Ranking 2 (P, Q,L\,L\, f(x,y) = y,y i+ i) 

-Ranking 2 (P, Q, L\,L % 2 , f(x, y) = y, yj. 

5 Si <— + Si-i. 

6 Find the smallest index w such that S w > k or w = m — 1 . 

7 r «- Ranking 2 (P, Q, Ljf, L™, /(x, y) = y, y TO ). 

8 return Selection 2 (P, Q, Lf, L%, f(x, y) = y, r - 1 + k - S w -i). 



Figure 10: The algorithm for the Minkowski Sum Selection problem with A > 2 linear constraints 
an a linear objective function. 
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